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Abstract. Nuclear parton distribution functions (NPDFs) are determined by a global analysis of experimental measurements 
on structure-function ratios F2 and Drell-Yan cross section ratios (y^y/^oY' Iheii uncertainties are estimated by the 
Hessian method. The NPDFs are obtained in both leading order (LO) and next-to-leading order (NLO) of a,. As a result, 
valence-quark distributions are relatively well determined, whereas antiquark distributions at x > 0.2 and gluon distributions 
in the whole x region have large uncertainties. The NLO uncertainties are slightly smaller than the LO ones; however, such a 
NLO improvement is not as significant as the nucleonic case. 
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INTRODUCTION 

From measurements on structure-function ratios /F2, 
we found that nuclear parton distribution functions 
(NPDFs) are not equal to corresponding nucleonic PDFs. 
Although there are many analyses on the PDFs in the nu- 
cleon, the NPDFs have not been investigated extensively. 
However, such studies are becoming more and more im- 
portant in recent years in order to precisely understand 
measurements of heavy-ion reactions at RHIC and LHC. 
It should lead to a better understanding on properties of 
quark-hadron matters. 

Such investigations are also valuable for applications 
to neutrino reactions. For an accurate determination of 
neutrino oscillation, nuclear corrections in the oxygen 
nucleus are important. Although the oscillation experi- 
ments are done in a relatively low-energy region, gross 
properties of cross sections could be described by using 
quark-hadron duality. In future, nuclear structure func- 
tions will be investigated by the MINER vA project and 
at neutrino factories. 

The unpolarized PDFs in the nucleon have been inves- 
tigated for a long time, and they are relatively well deter- 
mined. Uncertainties of the determined PDFs were also 
estimated recently, and such studies were extended to the 
polarized PDFs [ 1 ] and fragmentation functions [2] . On 
the other hand, we have been investigating parametriza- 
tion of the NPDFs and their uncertainties in a simi- 
lar technique [3, 4] by analyzing experimental data on 
nuclear structure-function ratios Fj* /F2 and Drell-Yan 
cross section rations g^y/^^dy- However, the uncertainty 



estimation was limited to the leading order (LO) in the 
previous analysis [4]. It is the purpose of this work to 
show the NPDFs and their uncertainties in both LO and 
next-to-leading order (NLO) and to discuss NLO im- 
provements [5]. 

ANALYSIS METHOD 

First, the functional form of our NPDFs is introduced. 
Since nuclear modifications are generally within the 
10%— 30% range for medium and large size nuclei, it 
is easier to investigate the modifications than the abso- 
lute NPDFs. We define the NPDFs at the initial scale 
(= Ql) as 

ft{x,Ql)^wt{x,A,Z)f,{x,Ql), (1) 

where //*(^, g^) is the parton distribution with type / 
(= My, dy, u, d, s, g) in a nucleus, /,(x,2q) is the cor- 
responding parton distribution in the nucleon, A is mass 
number, and Z is the atomic number. The variable 
is given = —q^ with the virtual photon momentum 
q in lepton scattering, and the Bjorken variable x is de- 
fined hy x = Q^/ (2Mv) with the energy transfer v and 
the nucleon mass M. We call vv,(x,A,Z) a weight func- 
tion, which indicates a nuclear modification for the type-! 
distribution. The weight functions are expressed by 

f A 7\ 1 , 1 \ai + bix + cix^+dix'^ 



where a, a,, bi, Ci, du and /3, are parameters to be deter- 
mined by a analysis. Here, the valence up- and down- 
quark parameters are the same except for and a^/,,. 
Since there is no data to find flavor dependence in an- 
tiquark modifications, the weight functions of u, d, and 
s are assumed to be the same at Qq. We impose three 
conditions, baryon-number, charge, and momentum con- 
servations, so that three parameters are fixed. 

The initial scale of the NPDFs is chosen Qq = 1 GeV^. 
They are evolved to experimental points. From the 
threshold at ^ m^., charm-quark distributions appear 
due to the evolution. Using these NPDFs, we calcu- 
late F2 /F2 and Drell-Yan ratios CT^y/^dy- param- 
eters are determined so as to minimize the total 
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where Rj indicates F2 /F2 and CJ^y/ CJ^y. The uncertain- 
ties of the determined NPDF are estimated by the Hes- 
sian method: 

[^fi^)?-^Z^Y.'-^H.^'J^, (4) 



where Hij is the Hessian matrix, ^, is a parameter, and 

t, indicates the optimum parameter set. The A.x^ value 
is taken as 13.7 so that the confidence level P becomes 
the one-(7-error range (P = 0.6826) for thirteen param- 
eters by assuming the normal distribution in the multi- 
parameter space. 



RESULTS 

We used data with > IGeV^. They consist of 290, 
606, 293, and 52 data points for F^/F{, F^/F^, 

F2/F2 {A' ^ D), and CToy/^^dy^ respectively. For the 
total 1241 data, we obtained the minimum values, 
xliJd.o.f.=1.35 and 1.21 for the LO and NLO. 

As examples, we show actual data with the theoretical 
LO and NLO ratios for the calcium nucleus together 
with their uncertainties in Figs. 1 and 2. The theoretical 
curves and their uncertainties are calculated at — 
10 GeV2 and = 50 GeV^ for the F2 and the Drell-Yan, 
respectively; however, the experimental data were taken 
at various values. There are discrepancies from the 
experimental data at x < 0.01 in Fig. 1, but they should 
be attributed to the difference [5]. The figures indicate 
good agreement between the data and the theoretical 
curves, and most data are within the uncertainties. As 
expected, the NLO uncertainties are smaller than the LO 
ones, especially at small x, but they are similar in the 
region, x > 0.02. As for Drell-Yan ratios in the region 
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FIGURE 1. Experimental data of Fj'^/F^ are compared 
with theoretical ratios, which are calculated at 2^ = 10 GeV^. 
The dashed and solid curves indicate LO and NLO results, 
the uncertainties are shown by the shaded bands. 
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FIGURE 2. Experimental data of a'^"/G^ are compared 
with the theoretical ratios calculated at = 50 GeV". The 
dashed and solid curves indicate LO and NLO results, and the 
uncertainties are shown by the shaded bands. 



X > 0.04, the NLO effects on the uncertainties are not 
obvious. 

In order to illustrate the nuclear dependence, we show 
the weight functions for all the analyzed nuclei and '^O 
at 2^ = 1 GeV^ in Fig. 3. The NPDFs in the oxygen nu- 
cleus are shown because they are important for neutrino- 
oscillation studies. As the mass number becomes larger 
in the order of D, '^He, Li, • ■ ■ , and Pb, the curves deviate 
from the line of unity. 

We provide a code for calculating the NPDFs and 
their uncertainties at our web site [6]. By supplying the 
kinematical conditions, x and Q^, and a nuclear species, 
one can obtain the NPDFs (m^, d"^, /, O^, d^, s^, c^, and 
^) numerically. The technical details on its usage are 
explained in Refs. [4, 5] and within the subroutine. 

Nuclear modifications and their uncertainties are 
shown for the calcium nucleus in Fig. 4. The valence- 
quark distributions are well determined in the medium 
and large x regions because the F^ structure functions 
are dominated by the valence-quark distributions and 
because the F2 ratios are accurately measured. At small 
X, the valence-quark modifications have also small 
uncertainties because of the baryon-number and charge 
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FIGURE 3. Nuclear modifications w,- (? = My, d,,, q, and g) 
are shown in the NLO for all the analyzed nuclei and '^O at 
= 1 GeV2. 



conservations. 

The antiquark distributions are well determined at x < 
0.2 due to the F2 and Drell-Yan data. However, they have 
large uncertainties at jc > 0.2 because there is no Drell- 
Yan data which constrains the antiquark modifications. 
Future Drell-Yan measurements at large x should im- 
prove the situation of the antiquark distributions. There 
are J-PARC, Fermilab-E906, and GSI-FAIR projects, in 
which Drell-Yan processes will be investigated. 

We found that the uncertainty bands for the gluon are 
very large, which indicates that the gluon modifications 




3 0.8 Q^=lGeV^ 
^ ^ LO 



0.6 
0.4 



NLO 




0.001 0.01 0.1 1 
JC 

FIGURE 4. Nuclear modifications of the PDFs and their un- 
certainties are shown for the calcium nucleus at 2^ = 1 GeV^. 
The dashed and solid curves indicate LO and NLO results, and 
their uncertainties are shown by the shaded bands. 



cannot be well determined in the whole x region. The 
gluon distributions contribute to the F2 and Drell-Yan 
ratios as higher-order effects. Therefore, the gluon dis- 
tributions cannot be accurately determined especially in 
the LO analysis. Some improvements are expected in the 
NLO analysis. In fact. Fig. 4 indicates that the NLO un- 
certainty band for the gluon becomes smaller than the 
LO one. However, it is not as clear as the improvements 
in the polarized PDFs [1] and the fragmentation func- 
tions [2]. The gluon distribution in the nucleon has been 
determined mainly by using dependence of the struc- 
ture function F2. In the nuclear case, the dependen- 
cies of the ratios F2 /F2 are not accurately measured. It 
leads to the large uncertainty bands in the gluon modifi- 
cations even in the NLO analysis. We hope that accurate 
data will be provided at future electron facilities such as 
eRHIC and eLIC. 

Nuclear modifications in the deuteron are also inves- 
tigated in our recent analysis [5]. Since the deuteron 
data are used for determining the "nucleonic" PDFs after 
some nuclear corrections, current PDFs in the nucleon 
could partially contain nuclear effects. Proper nuclear 
corrections should be taken into account in the PDF anal- 
ysis of the nucleon to exclude such effects. 



SUMMARY 

By the global analyses of the data on the nuclear F2 
and Drell-Yan ratios, the NPDFs have been determined 
in both LO and NLO. Although the valence-quark dis- 
tributions and antiquark ones at x < 0.2 are well deter- 
mined, the uncertainty bands are very large in the an- 
tiquark modifications at jc > 0.2 and the gluon ones at 
whole X even in the NLO analysis. We need future exper- 
imental efforts for determining all the nuclear modifica- 
tions. Our NPDFs and their uncertainties can be calcu- 
lated by the code supplied in Ref. [6]. 
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